Interplay of halogen bonding and solvation in protein-ligand binding.

Halogen bonding is increasingly utilized in efforts to achieve high affinity and selectivity of molecules designed to bind proteins, making it paramount to understand the relationship between structure, dynamics, and thermodynamic driving forces. We present a detailed analysis addressing this problem using a series of protein-ligand complexes involving single halogen substitutions - F, Cl, Br, and I - and nearly identical structures. Isothermal titration calorimetry reveals an increasingly favorable binding enthalpy from F to I that correlates with the halogen size and σ-hole electropositive character, but is partially counteracted by unfavorable entropy, which is constant from F to Cl and Br, but worse for I. Consequently, the binding free energy is roughly equal for Cl, Br, and I. QM and solvation-free-energy calculations reflect an intricate balance between halogen bonding, hydrogen bonds, and solvation. These advances have the potential to aid future drug design initiatives involving halogenated compounds.

Cosolvent Dimethyl Sulfoxide Influences Protein-Ligand Binding Kinetics via Solvent Viscosity Effects: Revealing the Success Rate of Complex Formation Following Diffusive Protein-Ligand Encounter.

Protein-ligand-exchange kinetics determines the duration of biochemical signals and consequently plays an important role in drug design. Binding studies commonly require solubilization of designed ligands in solvents such as dimethyl sulfoxide (DMSO), resulting in residual amounts of DMSO following titration of solubilized ligands into aqueous protein samples. Therefore, it is critical to establish whether DMSO influences protein-ligand binding. Here, we address the general and indirect effect of DMSO on protein-ligand binding caused by solvent viscosity, which is strongly dependent on the relative concentrations of DMSO and water. As a model system, we studied the binding of a drug-like ligand to the carbohydrate recognition domain of galectin-3 in the presence of variable amounts of DMSO. We used isothermal titration calorimetry to characterize binding thermodynamics and 15N NMR relaxation to monitor kinetics. The binding enthalpy is not affected, but we observe a subtle trend of increasingly unfavorable entropy of binding, and consequently decreased affinity, with increasing DMSO concentration. The increasing concentration of DMSO results in a reduced association rate of binding, while the dissociation rate is less affected. The observed association rate is inversely proportional to the viscosity of the DMSO-water mixture, as expected from theory, but significantly reduced from the diffusion-controlled limit. By comparing the viscosity dependence of the observed association rate with that of the theoretical diffusion-controlled association rate, we estimate the success rate of productive complex formation following an initial encounter of proteins and ligands, showing that only one out of several hundred binding "attempts" are successful.

Bromodomain Interactions with Acetylated Histone 4 Peptides in the BRD4 Tandem Domain: Effects on Domain Dynamics and Internal Flexibility.

The bromodomain and extra-terminal (BET) protein BRD4 regulates gene expression via recruitment of transcriptional regulatory complexes to acetylated chromatin. Like other BET proteins, BRD4 contains two bromodomains, BD1 and BD2, that can interact cooperatively with target proteins and designed ligands, with important implications for drug discovery. Here, we used nuclear magnetic resonance (NMR) spectroscopy to study the dynamics and interactions of the isolated bromodomains, as well as the tandem construct including both domains and the intervening linker, and investigated the effects of binding a tetra-acetylated peptide corresponding to the tail of histone 4. The peptide affinity is lower for both domains in the tandem construct than for the isolated domains. Using 15N spin relaxation, we determined the global rotational correlation times and residue-specific order parameters for BD1 and BD2. Isolated BD1 is monomeric in the apo state but apparently dimerizes upon binding the tetra-acetylated peptide. Isolated BD2 partially dimerizes in both the apo and peptide-bound states. The backbone order parameters reveal marked differences between BD1 and BD2, primarily in the acetyl-lysine binding site where the ZA loop is more flexible in BD2. Peptide binding reduces the order parameters of the ZA loop in BD1 and the ZA and BC loops in BD2. The AB loop, located distally from the binding site, shows variable dynamics that reflect the different dimerization propensities of the domains. These results provide a basis for understanding target recognition by BRD4.

Entropy-Entropy Compensation between the Protein, Ligand, and Solvent Degrees of Freedom Fine-Tunes Affinity in Ligand Binding to Galectin-3C.

Molecular recognition is fundamental to biological signaling. A central question is how individual interactions between molecular moieties affect the thermodynamics of ligand binding to proteins and how these effects might propagate beyond the immediate neighborhood of the binding site. Here, we investigate this question by introducing minor changes in ligand structure and characterizing the effects of these on ligand affinity to the carbohydrate recognition domain of galectin-3, using a combination of isothermal titration calorimetry, X-ray crystallography, NMR relaxation, and computational approaches including molecular dynamics (MD) simulations and grid inhomogeneous solvation theory (GIST). We studied a congeneric series of ligands with a fluorophenyl-triazole moiety, where the fluorine substituent varies between the ortho, meta, and para positions (denoted O, M, and P). The M and P ligands have similar affinities, whereas the O ligand has 3-fold lower affinity, reflecting differences in binding enthalpy and entropy. The results reveal surprising differences in conformational and solvation entropy among the three complexes. NMR backbone order parameters show that the O-bound protein has reduced conformational entropy compared to the M and P complexes. By contrast, the bound ligand is more flexible in the O complex, as determined by 19F NMR relaxation, ensemble-refined X-ray diffraction data, and MD simulations. Furthermore, GIST calculations indicate that the O-bound complex has less unfavorable solvation entropy compared to the other two complexes. Thus, the results indicate compensatory effects from ligand conformational entropy and water entropy, on the one hand, and protein conformational entropy, on the other hand. Taken together, these different contributions amount to entropy-entropy compensation among the system components involved in ligand binding to a target protein.

Rapid measurement of heteronuclear transverse relaxation rates using non-uniformly sampled R1ρ accordion experiments.

Multidimensional, heteronuclear NMR relaxation methods are used extensively to characterize the dynamics of biological macromolecules. Acquisition of relaxation datasets on proteins typically requires significant measurement time, often several days. Accordion spectroscopy offers a powerful means to shorten relaxation rate measurements by encoding the "relaxation dimension" into the indirect evolution period in multidimensional experiments. Time savings can also be achieved by non-uniform sampling (NUS) of multidimensional NMR data, which is used increasingly to improve spectral resolution or increase sensitivity per unit time. However, NUS is not commonly implemented in relaxation experiments, because most reconstruction algorithms are inherently nonlinear, leading to problems when estimating signal intensities, relaxation rate constants and their error bounds. We have previously shown how to avoid these shortcomings by combining accordion spectroscopy with NUS, followed by data reconstruction using sparse exponential mode analysis, thereby achieving a dramatic decrease in the total length of longitudinal relaxation experiments. Here, we present the corresponding transverse relaxation experiment, taking into account the special considerations required for its successful implementation in the framework of the accordion-NUS approach. We attain the highest possible precision in the relaxation rate constants by optimizing the NUS scheme with respect to the Cramér-Rao lower bound of the variance of the estimated parameter, given the total number of sampling points and the spectrum-specific signal characteristics. The resulting accordion-NUS R1ρ relaxation experiment achieves comparable precision in the parameter estimates compared to conventional CPMG (Carr-Purcell-Meiboom-Gill) R2 or spin-lock R1ρ experiments while saving an order of magnitude in experiment time.

Backbone 1H, 13C, and 15N resonance assignments of BoMan26A, a ß-mannanase of the glycoside hydrolase family 26 from the human gut bacterium Bacteroides ovatus.

Bacteroides ovatus is a member of the human gut microbiota. The importance of this microbial consortium involves the degradation of complex dietary glycans mainly conferred by glycoside hydrolases. In this study we focus on one such catabolic glycoside hydrolase from B. ovatus. The enzyme, termed BoMan26A, is a ß-mannanase that takes part in the hydrolytic degradation of galactomannans. The crystal structure of BoMan26A has previously been determined to reveal a TIM-barrel like fold, but the relation between the protein structure and the mode of substrate processing has not yet been studied. Here we report residue-specific assignments for 95% of the 344 backbone amides of BoMan26A. The assignments form the basis for future studies of the relationship between substrate interactions and protein dynamics. In particular, the potential role of loops adjacent to glycan binding sites is of interest for such studies.

Structural model of dodecameric heat-shock protein Hsp21: Flexible N-terminal arms interact with client proteins while C-terminal tails maintain the dodecamer and chaperone activity.

Small heat-shock proteins (sHsps) prevent aggregation of thermosensitive client proteins in a first line of defense against cellular stress. The mechanisms by which they perform this function have been hard to define due to limited structural information; currently, there is only one high-resolution structure of a plant sHsp published, that of the cytosolic Hsp16.9. We took interest in Hsp21, a chloroplast-localized sHsp crucial for plant stress resistance, which has even longer N-terminal arms than Hsp16.9, with a functionally important and conserved methionine-rich motif. To provide a framework for investigating structure-function relationships of Hsp21 and understanding these sequence variations, we developed a structural model of Hsp21 based on homology modeling, cryo-EM, cross-linking mass spectrometry, NMR, and small-angle X-ray scattering. Our data suggest a dodecameric arrangement of two trimer-of-dimer discs stabilized by the C-terminal tails, possibly through tail-to-tail interactions between the discs, mediated through extended IXVXI motifs. Our model further suggests that six N-terminal arms are located on the outside of the dodecamer, accessible for interaction with client proteins, and distinct from previous undefined or inwardly facing arms. To test the importance of the IXVXI motif, we created the point mutant V181A, which, as expected, disrupts the Hsp21 dodecamer and decreases chaperone activity. Finally, our data emphasize that sHsp chaperone efficiency depends on oligomerization and that client interactions can occur both with and without oligomer dissociation. These results provide a generalizable workflow to explore sHsps, expand our understanding of sHsp structural motifs, and provide a testable Hsp21 structure model to inform future investigations.